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Abstract. In this paper, we consider dynamical behavior of astrophysical objects such as 
galaxies and dwarf galaxies taking into account both the gravitational attraction between 
them and the cosmological expansion of the Universe. First, we obtain the general system 
of equations and apply them to some abstract systems of galaxies. Then we investigate 
the collision between the Milky Way and Andromeda in future. We demonstrate that for 
currently known parameters of this system the collision is hardly possible because of the 
angular momentum. These galaxies will approach the minimum distance of about 290 Kpc 
in 4.44 Gyr from present, and then begin to run away irreversibly from each other. We also 
define the region in the vicinity of our Local Group where the formation of the Hubble flows 
starts. For such processes, the zero-acceleration surface (where the gravitational attraction 
is balanced by the cosmological accelerated expansion) plays the crucial role. We show that 
such surface is absent for the Local Group. Instead, we find two points and one circle with 
zero acceleration. Nevertheless, there is a nearly closed area around the MW and M31 where 
the absolute value of the acceleration is approximately equal to zero. The Hubble flows are 
formed outside of this area. 
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1 Introduction 

The progress in modern observational cosmology at scales much smaller than the cell of 
uniformity size (see, e.g., [1-6]) enables to use the new observational data to test different 
cosmological models. With the help of these data, we can reconstruct the history of galaxies, 
groups and clusters of galaxies as well as to predict their future. For example, we can explain 
the formation of the Hubble flows in the vicinity of the group of galaxies [7] or predict a 
possible collision of the Milky Way and Andromeda in future [8] . 

According to recent astronomical observations, there is no clear evidence of spatial 
homogeneity up to sizes ~ 150 Mpc [9]. Deep inside of such scales and on late stages of 
evolution, the Universe consists of a set of discrete inhomogeneities (galaxies, groups and 
clusters of galaxies) which disturb the background Friedmann Universe. Hence, classical 
mechanics of discrete objects provides more adequate approach than hydrodynamics with 
its continuous flows. In our previous paper [10], we have elaborated this approach for an 
arbitrary number of randomly distributed inhomogeneities on the cosmological background 
and found a gravitational potential of this system. We have shown that this potential has 
the most natural form in the case of the Friedmann-Robertson- Walker metrics with the 
hyperbolic space. Therefore, having the gravitational potential of an arbitrary system of 
inhomogeneities, we can investigate their motion taking into account both the gravitational 
attraction between them and the cosmological expansion of the Universe. 

In the present paper, we continue this investigation. First, we obtain the general system 
of equations of motion for such system and apply these equations to abstract groups of 
galaxies to show the effects of gravitational attraction and cosmological expansion. Then, we 
consider our Local Group to investigate the mutual motion of the Milky Way and Andromeda. 
Contrary to the conclusions of the paper [8], we show that for currently known parameters 
of this system, the collision is hardly possible in future because of the angular momentum. 
These galaxies will approach the minimum distance of about 290 Kpc in 4.44 Gyr from 
present, and then begin to run away irreversibly from each other. Additionally, we define the 
region in the vicinity of our Local Group where the formation of the Hubble flows starts. For 
such processes, the zero-acceleration surface (where the gravitational attraction is balanced 



- 1 - 



by the cosmological accelerated expansion) plays the crucial role. We take into account the 
geometry of the system consisting of two giant galaxies (MW and M31) at the distance 0.78 
Mpc. Obviously, if this surface exists, it does not have a spherical shape. We show that 
such surface is absent for the Local Group. Instead, we find two points and one circle with 
zero acceleration. Nevertheless, there is a nearly closed area around the MW and M31 where 
the absolute value of the acceleration is approximately equal to zero. The Hubble flows are 
formed outside of this area. 

The paper is organized as follows. In section 2 we obtain the general system of equations 
of motion for arbitrary distributed inhomogeneities in the open Universe. We apply these 
equations to abstract systems of galaxies consisting of three and four galaxies in section 3. 
In section 4, we investigate the mutual motion of the Milky Way and Andromeda in future. 
In section 5, we define the zero-acceleration region for our Local Group. The main results 
are summarized in concluding section 6. 



2 General setup 

In our recent paper [10], we have shown that the "comoving" gravitational potential for a 
system of gravitating masses m; is 

^ exp(-2Zj) AttGnp 

i 

where Gn is the Newtonian gravitational constant, 7i = const is the comoving average rest 
mass density and Zj denotes the comoving geodesic distance between the i-th mass rrii and 
the point of observation in the open Universe, i.e. in the hyperbolic space. This formula has 
a number of advantages with respect to the flat and spherical space cases [10]. First, this 
potential is finite at any point of space (excluding, of course, the positions of the particles with 
li = 0). Second, the presence of the exponential function enables us to avoid the gravitational 
paradox (the Neumann-Seeliger paradox). The 7>term does not spoil this property because 
the averaged gravitational potential Tp is equal to zero. Third, the gravitating masses can 
be distributed completely arbitrarily. It is worth noting that, from different reasonings, the 
arguments in favour of the open Universe were also provided in the recent paper [11]. 

We consider the potential (2.1) against the cosmological background. In the ACDM 
model, the scale factor of the Universe reads [10] 



\lu_ 



2/3 

, (2.2) 



where we introduced dimensionless variables a = a/ao, i = H^t, and ao, Hq are the values 
of the scale factor a and the Hubble "constant" H = a/a = (da/dt)/a at the present time 
t = to (without loss of generality, we can put to = 0) . The standard density parameters are 



Ac 2 



where k = 8itGn/c 4 and A is the cosmological constant. It is worth noting that the solution 
(2.2) is common for any value of the curvature parameter /C because the density parameter for 
the curvature \Qjc\ = l^-l c2 /( a o-^o) ^ 1 l^] (where K = —1,0,-1-1 for open, flat and closed 
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Universes, respectively), and we can drop the spatial curvature term from the Friedmann 
equation. According to the seven-year WMAP observations [12], Hq ~ 70km/sec/Mpc ~ 
2.3 x lO-^sec" 1 ps (13.7 x lO 9 ) -1 ^ -1 , fl M « 0.27 and n A » 0.73. 

Below, we shall consider astrophysical objects (galaxies and their groups) deep inside 
of the cell of uniformity, i.e. for physical distances R < 150 Mpc. As we have shown in [10], 
the comoving distances in the cell of uniformity are much less than 1: Zj <C 1. For such small 
distances, we can use the Cartesian coordinates. Then, eq. (2.1) reads 

^-^ r — r j 3 

i 

where is the comoving radius- vector of the i-th gravitating mass rrij. Its Lagrange function 
is [10] 

m,<A' mja 2 v} 

d = — — + -hr 1 - > ( 2 - 5 ) 

where 

is the gravitational potential created by all the remaining masses at the point r = rj. In eq. 
(2-5), 

v\ = x? + y 2 + z 2 , (2.7) 

and Vi is the comoving peculiar velocity of the i-th gravitating source. It has the dimension 
(time) -1 . We would remind that we work in the weak-field limit where physical peculiar 
velocities are much less than the speed of light: avi <C c. 

It can be easily verified that, with respect to the physical coordinates 

Xi = axi, Yi = am, Zi = az t , (2.8) 

the Lagrange function is 



Ll ~ a + 2a? 



X;a — aXi ) + ( Y;a — aY; ) + ( Z,-,a — aZ, 



(2.9) 



where 

-d Ain. > 

(2.10) 

and the Lagrange equations for the i-th mass take the form 



<p i = -G N a 1 £ i m _ 3 \Ki - Kj\ = ^{X, i - Xj) 2 + (Y - Y 3 f + (Z { - Zjf , 



G N ^ m ^ Xi -y = i (X ia - aX) , (2.11) 



mjjXj-Xj) 1 

[(Xi - Xj) 2 + (Yi - Y 3 f + (Zi - Zjj- 

-G N £ 7 ^OTi-Yj) = l r fia _ &Y \ (212) 

i# [(X % - Xj) 2 + (Yi - Y 3 ) 2 + (Zi - Zj) 2 \ 

-G N £ = ^ ^ ~ **) • (2 - 13) 

& [(Xi - Xj) 2 + (Yi - Y 3 ) 2 + (Zi - Zj) 2 \ 

Now, we can apply these equations to real astrophysical systems such as a group or cluster 
of galaxies. To illustrate this, we consider first a number of abstract simplified examples. 
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3 Illustrative examples 



In this section, we consider simplified examples where all gravitating masses are on the same 
plane Z = 0, i.e. all Zi = 0. Then, eqs. (2.11)-(2.13) take the form 

d 2 Xi 1 v-^ mdXi-Xi) ld 2 a~ , . 

+ ^-^X t , i,j = l,...,N, (3.1) 



dt 2 m 



j^.[(X i -X j ) 2 + (Y i -Y j ) 2 ]y 2 ~adt 2 



d 2 Yi 1 v-^ mAYi-YA ^d 2 a~ Ar , , 

-^r = > — = J - J - + ^Yi, i,j = l,...,N, (3.2) 

(it 2 [(Xi-Xtf + iYi-Yj)*]*/* adt 2 y J 



H 2 \ 1/3 _ Xi /1O 12 M ^ 1/3 



where is a total number of masses and we introduced the dimensionless variables 

Xi = Xi 

Y = Yi 



Gattti J 0.95Mpc \ m 
H 2 \ 1/3 _ Y t /1O 12 M 



G^m J 0.95Mpc \ m 
t = H * t= 13.7 x 10V 



t = H t= — * q (3.3) 



and m = X^i 771 */-^ ^ s ^he average mass of the system. Obviously, the first terms in the 
right hand side of eqs. (3.1) and (3.2) are due to the gravitational attraction between masses 
and the second terms originate from the cosmological expansion of the Universe, which is 
described by eq. (2.2). We consider the stage of the accelerated expansion, i.e. d 2 d/di 2 > 0. 
Therefore, the competition between these two mechanisms defines the dynamical behavior 
of the masses. Depending on the initial conditions, they either collide with or move off each 
other. Let us demonstrate this with two particular examples. 

3.1 Three gravitating masses: N = 3 

Here, we study dynamics of three gravitating masses (N = 3) for two different cases. First, 
we consider the case with the initial coordinates (0,-1.5), (1.5,0) and (—1.5,1.5) for the 
first (i = 1), second (i = 2) and third (i = 3) gravitating masses, respectively, and with 
zero initial velocities 1 dXi/dt\^ =0 = 0, dYi/di\i =a = 0, i = 1,2,3. For simplicity, we also take 
mi = rri2 = = m. The numerical calculation of eqs. (3.1) and (3.2) shows (see figure 
1, the left panel) that in this case the cosmological expansion prevails over the gravitational 
attraction for all three masses and all of them move off each other. The red (inner) solid 
triangle corresponds to the initial positions at the initial moment t = 0. Then, we depict 
green (middle) and blue (outer) triangles at the moments t = 1 and t = 2, respectively. 
Solid triangles take into account both cosmological expansion and gravitational attraction, 
while the dashed triangles correspond to pure cosmological expansion. Obviously, the dashed 
triangles are similar but the gravitational attraction spoils this similarity (see the solid ones). 
Additionally, the solid triangles are inside of the corresponding dashed ones because the 
gravitational attraction slows the recession. 

In the second case, the initial coordinates are (0,-1), (1,0) and (—1,1), respectively, 
with the same zero initial velocities at t = 0. Then, the results of the numerical calculations 



1 The components of the velocity are given by the formula Xi = a±i + HXi. Therefore, to get zero value 
for this expression at some moment to, the peculiar velocity should compensate the Hubble velocity at this 
moment. 
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for times t = 1 (green triangles) and t = 2 (blue triangles) are shown in figure 1, the right 
panel. Similar to the previous case, solid (dashed) triangles correspond to the presence 
(absence) of the gravitational attraction. Therefore, the dashed triangles form the similar 
triangles. The solid green and blue triangles demonstrate that the gravitational attraction 
between first (i = 1) and second (i = 2) masses prevail the cosmological expansion, and these 
masses approach each other. However, the cosmological expansion dominates for the third 
(i = 3) mass and this mass moves away from the other two. 




Figure 1. Dynamics of three gravitating masses with zero initial velocities. The solid red triangles 
describe the initial positions at t = 0. The green and blue triangles correspond to the positions 
at the moments t = 1 and t — 2, respectively. Solid green and blue triangles take into account 
both cosmological expansion and gravitational attraction, while the corresponding dashed triangles 
disregard this attraction. Depending on the initial conditions, the gravitating masses move away from 
each other because the cosmological expansion prevails the gravitational attraction (the left panel), 
or some of masses can collide with each other in the case of the prevalence of the attraction (the right 
panel) . 



3.2 Four gravitating masses: N = 4 

Here, we study dynamical behavior of four gravitating masses (N = 4). The results of numer- 
ical calculations of eqs. (3.1) and (3.2) are depicted in figure 2. The left panel demonstrates 
the situation when the cosmological expansions prevails over the gravitational attraction and 
all masses move away from each other. The initial coordinates here are (1.5, 1.5), (1.5, —1.5), 
(—1.5,-1.5) and (—1.5,1.5). The right panel corresponds to the opposite case when the 
gravitational attraction dominates and masses approach each other. Here, the initial coor- 
dinates are (0.5,0.5), (0.5,-0.5), (—0.5,-0.5) and (—0.5,0.5). The meaning of color and 
type of lines is the same as for the previous example, i.e. the red squares correspond to the 
initial positions of masses at t = and green and blue squares show their positions at times 
t = 1, t = 2 (the left panel) and t = 0.3, t = 0.6 (the right panel), respectively. Solid green 
and blue squares take into account both cosmological expansion and gravitational attraction, 
while the corresponding dashed squares disregard this attraction. In contrast to the previous 
three-mass example, here, we alow the masses to rotate clockwise setting the following initial 
velocities (dXi/dt,dYi/dt) i=0 : (0,-1), (-1,0), (0,1), (1,0) (the left panel) and (0,-0.25), 



- 5 - 



(—0.25,0), (0,0.25), (0.25,0) (the right panel), respectively. Additionally, the orange curved 
line depicts the trajectory of one of the masses. 



Figure 2. Dynamics of four gravitating masses with non-zero initial velocities. The solid red squares 
correspond to the initial positions of masses at t = 0, and green and blue squares show their positions 
at times t = 1, t = 2 (the left panel) and t = 0.3, t = 0.6 (the right panel). Solid green and blue squares 
take into account both cosmological expansion and gravitational attraction, while the corresponding 
dashed squares disregard this attraction. The orange curved line depicts the trajectory of one of the 
masses. The left panel demonstrates the situation when the cosmological expansion prevails over the 
gravitational attraction and all masses move away from each other. The right panel corresponds to 
the opposite case when the gravitational attraction dominates and masses approach each other. 



4 Collision between Milky Way and Andromeda 

Now, we want to apply our method to real astrophysical objects. For this purpose, we consider 
our local group of galaxies, which consists of two giant galaxies (our Milky Way (MW) and 
Andromeda (M31)) and approximately 40 dwarf galaxies. At the present time, these giant 
galaxies are located at the distance 0.78 Mpc and move towards each other with the speed 
120 km/sec [8]. Therefore, in future they may encounter. The collision time was estimated 
recently in the paper [8], where the authors used the hydrodynamic approach. They found 
that the average time for the first passage is 2.8 Gyr and for the final merger is 5.4 Gyr. It 
is of interest to estimate also this time using our mechanical approach. We consider these 
two galaxies as pointdike gravitating masses. Obviously, such approach is valid at distances 
greater than the sizes of galaxies. For these two galaxies, we can apply our method up to 
the separation distance of the order of 100 Kpc, where the process of merger starts [8] 2 . Due 
to the very small intergalactic/intragroup medium density (which is only 5-200 times the 
average density of the Universe [14, 15]) we also neglect it in our calculations. Moreover, the 
viscosity of this very rare medium is extremely small. Consequently, the galaxies move in 
this medium without any essential resistance. 



2 According to the simulations carried out in [13], this distance can be extended up to 120-150 Kpc. 
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From eq. (2.9), it can be easily seen that the two-particle Lagrange function for two 
gravitating masses/galaxies (marked as the points A and B) is 



m A m B , m A 
CAB = Gn \R a -K b \ + 2^ 



X A a - aX A ) + {Y A a - aYj^ + ( Z A a - aZ A 



+ 



+ 



m B 
2a 2 



X B a - aX B ) + {y B a - ixY B \ + ( Z B a - aZ B 



(4.1) 



Let us introduce the projections Lx,Ly and L z of the distance between these masses and 
the coordinates of the center of mass: 



X A - X B = Lx , 
m A X A + m B X B 

= An 



(4.2) 
(4.3) 



m A + m B 

and the similar expressions for Ly,L z and Yq,Zq. Therefore, the absolute value of the 



distance is \~R A — ~R B \ = \ L 2 X + L\ + L| = L > 0. Then, eq. (4.1) reads 



m A m B 1 
L-AB = <~rN ; h 



m A m B ( . 2 2 • . 2r 2 , 



L 2a 2 [ 771 + 771 £ 

+ Lya 2 - ILyLyaa + a 2 Ly + L 2 z a 2 - 2L z L z aa + a^L| ) + 
+ {m A + mg) 



X a - aX ) + ( Y a - oYq] + ( Z a - aZ 



(4.4) 



or, in spherical coordinates 
C-AB = Gn 



L 2 wia + m B \a 2 



m A m B , 1 m A m B (a # _ 2 a ^ + ^ + ^ 2 + ^ sin2 ^ 2 + 

a 



1 



+ ^2( TO A + r7i B ) 



X a - aX ) + (Yqcl - aY^] + ( Z a - aZ 



It can be easily verified that Xq satisfies the equation 

aXo — aXo = 

with the following solution: 

a 



Xq — HX + — ( A~n( in ) - -ffn^O(m) I ! 



(4.5) 



(4.6) 



(4.7) 



where an = a(tn),-ffn = H(to), X ^ = Xn(tn) and X ^ = A~n(tn) are values at the initial 
time t = in. Therefore, A*o satisfies asymptotically (with increasing a) the Hubble law. The 
same conclusion takes place for Yq and Zq. 

Let us investigate now the relative motion of the galaxies. For this motion, the Lagrange 
function is 



Cab = G N ^-^ + - ""7^ I - 2-LL + L z + i/> 
L 2 itla + rn B \a z 



m,Am B 1 m A m B fa 2 2 



2„;,2 



(4.8) 
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where without loss of generality we put = tt/2. Therefore, the Lagrange equation for the 
separation distance is 

- 2m M 2 a T ,. n . 

L = - GN U + ^ + a L ' (49) 
where we introduced the reduced mass, the average mass and the angular momentum: 

m A m B _ _ m A + m B 2 • _ 

/U, m = , llL ip — M = const . (4-10) 



tua + rriB 



The first term in the right hand side of (4.9) is due to the gravitational attraction, the second 
term is the centrifugal force and the third term originates from the cosmological expansion 
of the Universe. To integrate eq. (4.9), we rewrite it with respect to the dimensionless 
quantities similar to ones in (3.3): 



L = L 



H 2 \ 1/3 L /1O 12 M ^ 1/3 



G^m J 0.95Mpc \ m 

V = V[^-) 1/3 ^ V flL^) , (4.11) 
[ H G N mJ 67km/sec 1 = ' 



m 



m ( H n \ 1/3 l V ± /io 12 M ^ /3 



\i \G\rn? ) 0.95Mpc 67km/sec V, m 

where the radial velocity V = L —> V = dL/dt and the transverse velocity V± = Lip. Then, 
eq. (4.9) reads 

d 2 L 2 M 2 \d 2 ~a~ /A . 

W = -p + T? + l!P L - (4 ' 12) 

Now, we integrate eq. (4.12) for parameters corresponding to the galaxies MW and 
M31. The masses of MW and M31 are of the order of 1O 12 M and 1.6 x 10 12 M Q respectively 
[8] 3 . The separation distance at present time 4 t = to is Lq 0.78 Mpc — > Lq 0.753, and 
the galaxies approach each other with the radial velocity V$ ~ —120 km/sec — > Vq ~ —1.633 
[8]. 

First, let us consider briefly the case of the zero angular momentum M = 0. If a > 0, 
as it happens at the present stage of the Universe evolution, we can introduce a distance of 
zero acceleration L cr where L = 0: 

- - i i \ 1 ' 3 

L ^-{-=m) • < 4 - 13 > 

where the deceleration parameter q = —(l/H 2 )(d/a) = —{d 2 a/dt 2 )/a. In the ACDM model, 
we get from the Friedmann equations that at present time qo « — ~ —0.595, and 

the zero acceleration distance is L cr (to) ~ 1.5. If at t = to the relative velocity L(to) = 0, 
then gravitating masses run away from each other (collide with each other) in the future for 
L(to) > L cr (to) (L(to) < L cr (to)). Obviously, the separation distance Lo ~ 0.753 between 
MW and Andromeda is less than L cr (to). Additionally, they have the non-zero radial velocity 
towards each other. Therefore, they will collide with each other. The result of numerical 



3 We take these values because we want to compare our results with the conclusions of this paper. More 
recent publications [16, 17] indicate a little bit higher values of the mass of MW. 
4 Without loss of generality we may put to = 0, to = 0. 
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calculation of this collision for M = is shown in figure 3, the left panel. The solid blue line 
takes into account both gravitational attraction and cosmological expansion while the red 
dashed line disregards the cosmological expansion. This picture demonstrates that the effect 
of the expansion is very small for relative motion of MW and M31. For example, the time 
to collision (from present) is t w 0.2670 — > 3.68Gyr and t w 0.2636 — > 3.63Gyr for blue and 
red lines, respectively. We remind that our approach works up to the separation distance 
L~0.1-> 100 Kpc when the stage of the galaxy merger starts. 




Figure 3. These figures show the change with time of the separation distance between the Milky Way 
and Andromeda starting from the present (t — 0). The initial (i.e. at present time) separation distance 
and radial relative velocity are 0.78 Mpc and —120 km/sec, respectively (0.753 and —1.633 in dimcn- 
sionless units). The solid blue lines take into account both gravitational attraction and cosmological 
expansion, while the red dashed lines disregard the cosmological expansion. The transverse velocity 
is absent in the left panel. Here, the collision between galaxies takes palace in t w 0.267 3.68 
Gyr from present (the blue line). In the right panel, the transverse velocity is equal to 100 km/sec. 
Here, the collision is absent and the smallest separation distance is L ps 0.28 — > L pa 290 Kpc at time 
t ps 0.324 — > t « 4.44 Gyr from present. For both pictures, the effect of the cosmological expansion is 
very small for the considered period of time. 

Let us turn now to the case of the non-zero angular momentum. The observations 
indicate the proper motion of Andromeda perpendicular to our line of sight. This transverse 
velocity V±q = V±(to) is less than 200 km/sec [18]. In [19], the authors found an even smaller 
estimate: V±o ~ 100 km/sec. In our calculations, we will adhere to this value, and with 
the help of eq. (4.11) we get M ~ 1.029. It can be easily seen from eq. (4.12) that fall to 
the center is absent because of the centrifugal barrier (M ^ 0). Therefore, the collision of 
the galaxies is possible if the smallest separation distance between them (which corresponds 
to the turning point) is less than the merger distance 100 — 150 Mpc. The result of the 
numerical integration of eq. (4.12) is shown in figure 3, the right panel. The solid blue 
line takes into account both gravitational attraction and cosmological expansion, while the 
red dashed line disregards the cosmological expansion. Similar to the previous case, the 
effect of the cosmological expansion is very small for the considered period of time. This 
picture demonstrates that for the given transverse velocity, the smallest separation distance 
is L pa 0.28 — > 290 Kpc at time t pa 0.324 — > AAA Gyr from present. This distance is much 
bigger than the merger distance. Therefore, for the chosen initial conditions, the collision 
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between the Milky Way and Andromeda is absent. The collision may take place for a smaller 
transverse velocity. For example, if Vj_o ~ 60 km/sec, then the smallest separation distance 
is 100 Kpc, which can be sufficient to start the merger. 



5 Formation of Hubble flows in the vicinity of the Local Group 

To study the formation of the Hubble flows in the vicinity of our group of galaxies, we 
need to determine the spatial distribution of vectors of acceleration of astrophysical objects 
(e.g., dwarf galaxies) in the gravitational field of two giant galaxies taking into account 
the cosmological expansion of the Universe. Obviously, near the galaxies, the vector must be 
oriented in the direction of galaxies due to the gravitational attraction, and with the distance 
from galaxies he has to turn in the opposite direction due to the cosmological accelerated 
recession. 

Let us investigate this effect for our Local Group. For this purpose, we consider a test 
particle/dwarf galaxy in the gravitational field of Andromeda and Milky Way. We study the 
picture at present time when the separation distance between M31 and MW is L$ = 0.78 
Mpc, and we do not take into account the relative motion of these galaxies. Of course, we 
can also consider dynamical evolution of this system but this effect is out of the scope of this 
section. It can be easily seen from the Lagrange function (2.9) that the Lagrange equation 
for a test particle is 

1( V -V)=-^ + 4r--V (5.1) 
dt\ a j adK a 2 a 

or, equivalently, 



V-"R = -I* (5.2) 
a a aK 

If the gravitational field is absent (i.e. ip = 0) then this equation has the following solution: 

a„ const 

V = -R + , 5.3 

a a 

where the first term is the Hubble velocity and the second term is the peculiar velocity. 
In dimensionless variables (see (3.3) and (4.11)) and without peculiar velocity this solution 
reads 

adt H y ' 

In the case of our Local Group, the gravitational potential is (see eq. (2.10)): 

m A m B 
ip = ip A + ip B , (p a = -a&N _ R | , ¥>B = -aG N — — - R j , (5.5) 

where, similar to the previous section, we mark MW and M31 by letters A and B, respectively. 
Here, we omitted the constant term ~ 7j because it does not contribute to equations of motion. 

For numerical calculation, we rewrite eq. (5.2) in dimensionless variables (3.3) and 
(4.11) as 

~ dV 1 d 2 d~ I dip , * 

W = ^ = ^-R 1, (5.6) 

dt adt 2 a OR 

where we introduced additionally the dimensionless potential 

<f>= _ 2/3 1- ( 5 - 7 ) 

a {H G N m) 1 
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It makes sense to rewrite (5.6) in components. For example, for the X-component we have 



with 



1 dip rriA 



a dX m 



~ _ dV x _ 1 d 2 a g 1 dip 
x ~ dt ~ a dt 2 a dX 



X-X A 



(5i 



+ 



m B 



J 2 +(Y-Y A y +(z-Za' 
X-X B 



3/2 



X - X B ) + (? - YbY + (z - Zb' 



3/2 



(5.9) 



and similar for the Y- and ^-components. The absolute value of the dimensionless accelera- 
tion is 

dV 



W 



dt 



1 



{H$G N my/ s 



dV 



dt 



(5.10) 



Eqs. (5.6) and (5.8) demonstrate how the cosmological expansion competes with the 
gravitational attraction. There is a characteristic region defined by the condition 



a dt 2 



1 d<p 
ddR 



(5.11) 



where the gravitational attraction is balanced by the cosmological expansion. The gravita- 
tional attraction is stronger inside of this region (for smaller distances) and the cosmological 
expansion prevails this attraction outside of this area (for larger distances) where the Hubble 
flows are formed. It takes place both for positive and negative values of the cosmological 
acceleration d. In the case of the accelerated expansion d > 0, we have |W| ~ and this 
characteristic area is called the region of zero acceleration. 

For our Local Group consisting of two giant galaxies MW and M31, we choose the origin 
of coordinates in the barycenter of these galaxies and X— axis along the line connecting MW 
andM31. Therefore, Xa = L^rns / {jnA^r^B)-, Xb = —LomA/(mA+mB), and Ya = Za = 0, 
Yb = Zb = 0. Additionally, due to the rotational symmetry around the X— axis, it is 
sufficient to consider the plane Z = 0. The 3-D picture can be easily reconstructed by rotation 
around this axis. Therefore, we investigate the distribution of the test body acceleration in 
the plane Z = 0. For the masses MW and M31, we take values from the previous section: 
m A ~ 1O 12 M and m B ~ 1.6 x 1O 12 M . 

On figure 4, we depict the absolute value of the acceleration (5.10) of the test body in 
the plane Z = 0. This modulus decreases from large values near the positions of MW and 
M31 to nearly zero (a red area around the peaks), and then it begins to increase again with 
the distance from the barycenter. The red area describes approximately the region of zero 
acceleration. Figure 5 depicts the vector field of the acceleration (5.6). It demonstrates the 
turn of these vectors from the directions towards the MW and M31 (in the vicinity of these 
galaxies) to outside (with distance from the galaxies). The central region (near MW and 
M31) is empty because we cut off the vectors with the magnitude |W| > 3. The yellow and 
green lines correspond to the conditions W x = and W y = 0, respectively. This figure shows 
that the vectors change their directions in the vicinity of the region |W| ~ 0. 
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Figure 4. This plot shows the absolute value of the acceleration of dwarf galaxies in the Local Group. 
The red area around the peaks corresponds approximately to the zero acceleration region. 
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Figure 5. This figure shows the vector field of the dwarf galaxy acceleration W. These vectors are 
directed towards the MW and M31 (the black points) in the vicinity of the galaxies and turn out with 
the distance from the galaxies. The yellow and green lines correspond to the conditions W x = and 
W y = 0, respectively. 

To define more exactly the structure of the zero-acceleration surface, we draw figure 6. 
The yellow and green lines correspond to the conditions W x = and W y = 0, respectively. 
Black points define the positions of the Milky Way (the right point) and Andromeda (the 
left point). Red points are defined by the condition W x = W y = — >■ |W| =0. The right 
panel takes into account both the gravitational attraction and the cosmological expansion, 
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while the left panel disregards the cosmological expansion. Obviously, in the case of only the 
gravitational attraction (the left panel), we have only one zero acceleration point between 
MW and M31 which is the analog of the Lagrange point L\. Much more reach picture 
happens in the presence of the cosmological accelerated expansion which competes with the 
gravitational attraction (the right panel). However, this panel shows that, strictly speaking, 
the zero acceleration surface is absent. Here, we have two additional red points on the X— axis 
and two vertical points. Clearly, due to the rotational symmetry, the latter two points are 
just the section of a zero acceleration circle by the plane Z = 0. Nevertheless, we can speak 
about the approximate zero acceleration surface because the elliptic-like green and yellow 
lines are very close to each other and they define the region where |W| ~ 0. We can see 
also that there are two regions where this surface has a discontinuity. It is clear from the 
rotational symmetry that these two regions belong to the round belt (they are the section of 
this belt by the plane Z = 0). Therefore, inside of the approximate zero acceleration surface 
the gravitational attraction is stronger than the cosmological expansion while outside of this 
surface the cosmological expansion prevails over the gravitational attraction. Obviously, 
the Hubble flows are formed in the latter region. Additionally, we can see that there is 
an asymmetry in directions along X and Y axes. The characteristic distances from the 
barycenter to the zero acceleration surface are \X\ ~ 1.6 — > \X\ ~ 1.65 Mpc and |Y| ~ 
1.45 \Y\ » 1.5 Mpc. 



1 ■ m 










Figure 6. Here, we depict the contour plot of the absolute value of the acceleration |W| in our 
Local Group (left and right black points are M31 and MW, respectively). The yellow and green lines 
correspond to the conditions W x = and W y = 0, respectively. Red points define the positions of the 
zero acceleration: W x = W y = — > \W\ = 0. The right panel takes into account both the gravitational 
attraction and the cosmological expansion, while the left panel disregards the cosmological expansion. 
The elliptic-like green curve together with the neighboring yellow curve defines the region where 
|W| w (the right panel). 
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6 Conclusion 



In this paper, we have considered the motion of astrophysical objects deep inside of the cell 
of uniformity where both the gravitational attraction between them and the cosmological 
expansion of the Universe play the role. To describe this, we obtained the general system of 
equations of motion for arbitrary distributed inhomogeneities in the open Universe. To show 
the competitive effects of gravitational attraction and cosmological recession, we considered 
two illustrative abstract examples where the systems of galaxies consist of three and four 
galaxies. 

Then, we investigated our Local Group which consists of two giant galaxies the Milky 
Way and Andromeda and approximately 40 dwarfs galaxies. According to recent observa- 
tions, these giant galaxies move towards each other with the relative velocity ~ 100 km/sec 
and may collide in future. Such process was investigated, e.g., in [8] where the authors 
used the hydrodynamic approach. They found that the average time for the first passage 
is 2.8 Gyr and for the final merger is 5.4 Gyr. Contrary to the conclusions of this paper, 
our mechanical approach has shown that for currently known parameters of this system, the 
collision is hardly possible in future because of the angular momentum. These galaxies will 
reach the minimum distance of about 290 Kpc in 4.44 Gyr from present, and then begin to 
run away irreversibly from each other. 

We have also defined the region in the vicinity of our Local Group where the Hubble 
flows start to form. For such processes, the zero acceleration surface (where the gravitational 
attraction is balanced by the cosmological accelerated expansion) plays the crucial role. We 
have taken into account that two giant galaxies MW and M31 are located at the distance of 
0.78 Mpc from each other. Obviously, if this surface exists, it does not have a spherical shape 
for given geometry. We have shown that such surface is absent for the Local Group. Instead, 
we have found two points and one circle with zero acceleration. Nevertheless, there is the 
nearly closed area around the MW and M31 where the absolute value of the acceleration is 
approximately equal to zero. The Hubble flows are formed outside of this area. 
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